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ABSTRACT 

We present an extension of the scale invariance formalism by'Heinz & Sunvaev'('2003") which 
includes the effects of radiative cooling on the X-ray synchrotron emission from core domi- 
nated jets and outflows. We derive the scaling relations between the radio luminosity L^, the 
synchrotron X-ray luminosity Lx, the mass M, and the accretion rate fn of an accreting black 
hole. We argue t hat the inclusion of sy nchrotron losses into the scaling relations is required by 
the data used in 'Merloni et al. (2003) to define the "fundamental plane" correlation between 
Lr, ix, and M in AGNs amd Galactic X-ray binaries. We then fit the new scaling relations 
to the "fundamental plane" correlation. This allows us to derive statistical constraints on the 
contribution from jet synchrotron emission to the X-rays from black holes, or, alternatively, on 
the particle acceleration scheme at work in these jets. We find that, in order for jet synchrotron 
X-rays to be consistent with the observed — ix — M correlation, the global particle spec- 

/Q H +0.6 \ 

trum must be both steep: f{j) oc 7 ^ ' -« s' and unbroken down to energies well below the 
radio synchrotron regime. Such a steep particle distribution would imply steeper X-ray spectra 
than observed in most sources contributing to the correlation with avalailable X-ray spectral 
index measurements. We suggest that, unless some assumptions of the scale invariance hy- 
pothesis are broken, another emission mechanism (likely either disk radiation and/or inverse 
Compton scattering in the jet) contributes significantly to the X-rays of a sizeable fraction of 
the sources. 

Key words: radiative mechanisms: non-thermal - galaxies: active - galaxies: nuclei - galax- 
ies: jets - x-rays: binaries - radio: continuum: general 



1 INTRODUCTION 

For many years, the standard paradigm for the X-ray emission in 
accreting black hole systems has bee n that the X-rays are pro- 
duced in the central accretion flow (e.g.. Shakura & Sunvae'vll973l : 
[shaciro et al. 1976; Naravan & Yi 1994), while radio emission was 
traditionally attributed to a radio jet. Only in the most powerful, 
beamed jet sources was X-ray synchrotron and inverse Compton 
emission from the jet invoked to ex plain the featureless powerlaw 
spectra of Blazars and BL-Lacs (e.g. lFossati et alll998h . 

In recent years, however, this paradigm has been called into 
question. Beginning with a re-evaluation of the emission processes 
in Galactic X-ray binaries (XRBs), an alternative view has been 
proposed to interpret a large part of the spectrum (including X-rays) 
of black holes at moderate to low accretion rates entirely as due to 
the em ission from the jet <Markoff et alJl200lL EoObI : Ipalcke et"ai] 
|2004 . 

It has been difficult to find unique signatures that would allow 
us to distinguish between these two alternatives. Broad band spec- 
tral modeling is the most promising approach, but so far it has not 
been shown conclusively that either of the two approaches fails. 



The strong correlation between the radio and the X-ray emission 
in X-ray binaries that has been found recently iCorbel et all2003l 
iGallo et all2003h can be interpreted in terms of both pictures. Sim- 
ilarly, the correlation between radio luminosity Lr, X-ray luminos- 
ity Lx, and black hole m ass M is co nsistent with both pictures 
jMerloni et al. 2003.:.Falcke et aj2004) . 

In order to inves t igate the relation between Lxr, Lx, and 
M, iHeinz & SunvaevI i2003h derived robust and very general 
scaling relations between these observa bles which are in d epen- 
dent of the underlying jet model (e.g. Bla ndford & Ostrikeilll978|: 
iFalcke & Biermaiinll995h . provided that jet physics is scale in- 
variant under changes in black hole mass M and accretion rate 
771 = AL/MEdd (where A^Edd is the Eddington accretion rate). 
The striking agreement with the subsequentl y found "fundamen- 
tal plane" relation of accreting black holes iMerloni et alJl2003l 
iFalcke et all2004 seems to confirm the scale invariance hypothe- 
sis^. 

^ While the influence of interactions between the jet and the environment 
must affect the radiative signature of jets on larger scales 1 Heinz & Sunyaej 
l2002t) . we note that the scaling relation predicted by IHeinz & Sunvae\l 
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The original treatment of iHeinz & Sunvae^ i2003h did, how- 
ever, leave out some critical physical ingredients required to de- 
scribe the X-ray synchrotron emission from jets, namely, radiative 
cooling of the particle distribution, which alters the spectral shape 
at higher energies. In this paper, we will present a n extension of 
the formalism developed in lHeinz & Sunvaevli2003l) that allows us 
to include radiative cooling (32|and i|Al> . We then proceed to in- 
vestigate the data of the "fundamental plane" relation for signs of 
radiative cooling and discuss limits we can put on the contribution 
of X-ray synchrotron radiation from jets to the spectra of AGNs 
and XRBs by fitting the derived scaling relations to the fundamen- 
tal plane correlation (S|3)- Section|4|summarizes the conclusions we 
draw from this analysis. 



2 SCALE INVARIANCE IN JETS WITH RADIATIVE 
LOSSES 

In iHeinz & Sunvae^ i2003 h we derived the scaling relations be- 
tween synchrotron radio luminosity Lr, black hole mass M, and 
accretion rate Ad. Because radiative losses introduce another scale 
into the problem, we chose to ignore the effect of synchrotron cool- 
ing on the particle distribution, given the scope of the paper. As we 
will show below, it is possible to incorporate these effects without 
loss of generality, because the scale introduced by cooling depends 
itself only on the mass and accretion rate of the black hole and does 
not enter into the dynamical properties of the jet plasma at all. This 
is because cooling affects primarily the highest energy particles, 
while most of the jet kinetic energy and thrust are carried by the 
low energy particles (in other words, jets are believed to be radia- 
tively inefficient). 

2.1 The scale invariance hypothesis for cooled jets 

In order to include the effects of radiative cooling into the scale 
invariance formalism (Heinz & Sunvaev 2003), we shall briefly re- 
view the aspects of the scale invariance ansatz relevant to this arti- 
cle. For a detailed discussion, we refer the interested reader to the 
original article. 

The scale invariance ansatz assumes that there is only one rel- 
evant scale in the problem of jet dynamics (at least in the cores 
of jets which we are focusing on here): the gravitational radius 
rg = GM/c? . It also assumes that the structure of all core jets^ 
from black holes is the same when set in relation to rg. Mathe- 
matically, the scale invariance assumption can be formulated as the 
requirement that any dynamical quantity / relevant to calculating 
the synchrotron emission from jets (such as the normalization of the 
powerlaw distribution of electrons, the magnetic field, or the parti- 
cle pressure) can be decomposed into two functions (j> f [M, rh) and 
ipfix)' such that 

}{M,m,r) = 4,}{M,m)^f{x) (D 

where the definition x ~ ''/'"g implies that tpf depends only on 
scale-free coordinates (x, y, z) /rg. 

The functions 0/ define the normalization of / at the base 
of the jet (as provided by the accretion disk), while the structure 
functions describe the variation of / along the jet, relative to 

1200 3 l) for optically thin synchi'otron emission agrees veiy well with the 
tentative mass and accretion rate d ependence of the o ptically thin radio lu- 
minosity of extended jets found by iSams et alj Il99d) . 



its value at the base of the jet. Any non-trivial dependence of / 
on the black hole spin a is assumed to be implicit in the functions 
i/i / and will not be discussed further. For consistency, we will also 
assume that the jet velocity /JF is independent of black hole mass 
and accretion rate, such that /?F = ^/3r(x)- It is still not known 
what the true jet plasma velocities are, especially close to the base 
of the jet, but future constraints will hopefully shed some light on 
this subject. If this assumption is not satisfied, the formalism can 
easily be altered to include such an M or m dependence as long as 
the four-velocity (3V can be written in the form of eq. Q. 

Given these assumptions, we can proceed to include radiative 
cooling into the treatment. We will assume that a fraction of the jet 
particles are accelerated at some arbitrary time to (measured in the 
jet frame) into a powerlaw such that 

/(7,to) =<^C7"'' (2) 

and then evolve through adiabatic and radiative cooling as they 
move downstream from the acceleration region. 

While it is not necessary to assume that this acceleration takes 
place in one particular location (it can be spatially extended, in 
which case the integral in eq. <A15> can be performed in a piece- 
wise fashion for every acceleration zone), it will be necessary to 
make the assumption that the spatial distribution of the accelera- 
tion (i.e., the various positions of the acceleration regions measured 
in units of gravitational radii rg) and the shape of the powerlaw 
are invariant under changes of M and m (in line with the assump- 
tion s necessary to deriv e the original scale invariance properties in 
IHeinz & Sunvae^^2003^ . This implies that the injection time in the 
jet frame to is proportional to to oc rg cx M. This is reasonable if 
the particle acceleration is associated with a hydrodynamic process, 
such as a recoUimation shock, since the underlying hydrodynamics 
is also assumed to be scale invariant. 

Furthermore, we will assume that the particle pressure is dom- 
inated by the low energy end of the particle distribution, which 
implies that the slope of the particle distribution produced by the 
acceleration process is steeper than p J> 2. Classical Fermi I 
shock acceleration in stro ng, non-relativistic shocks pr o duces spec- 
tra with index p ^ 2 telandford & Ostrikeil IT973: lBelllll978t 
Drurv et al. '1982^, while the spectral index produced in strong rel- 
ativistic shocks is som ewhat steeper, p ~ 2.3 (Gallant et al. 199^; 
lAchterberg et all200lh . 

These indices agree nicely with the observed optically thin ra- 
dio synchrotron spectra of extended jets with spectral indices of 
order 0.5 <; Or ^ 0.7 (where the radio spectral index is defined as 
Qr = — 9(log F-c)/ 9(log v)). For this reason, we will assume in the 
following that the fiducial spectral index of the uncooled electron 
distribution falls into the range 2 ^ p ^ 2.3. The fiducial uncooled 
synchrotron spectral index falls into the range 0.5 ^ a ^ 0.65. 
More recent efforts studying particle acceleration in jets are focus- 
ing on continu ous particle acce leration mechanisms like reconnec- 
tion (e.g. iLarrabee et aljEoOSi) . which might be necessary to ex- 
plain the lack of synchrotron cooling signatures in some extended 
AGN jets like 3C273 (Jester et al. 2001.) . 

2.2 Scaling relations 

With the basic set of assumptions in place, we proceed to derive the 
scaling relations for the X-ray synchrotron emission from jet cores 
including the effects of radiative cooling on the particle distribu- 
tion. The complete mathematical derivation used in the remainder 
of this paper is presented in the appendix. For the casual reader, we 
will present a sketch of the derivation in the following paragraphs. 
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Figure 1. Sketch of viable synchrotron spectra in different regimes for three 
different measured radio-to-X-ray ratios. Optically thick parts are shown as 
solid lines, optically thin, uncooled parts are shown as dashed lines, op- 
tically thin, radiatively cooled parts are shown as dotted lines. The grey 
dashed line through Lx3 shows an X-ray spectrum that is flatter than the 
global spectrum, which would require a concave spectrum and would be 
inconsistent with simple synchrotron X-rays. 



Synchrotron cooling of a population of electrons in the pres- 
ence of a source of fresh powerlaw particles (which is the appro- 
priate model for a continuous jet) will produce a cooling break in 
the global spectrum, as sketched in Fig.Q Given the field strength 
B, the characteri stic "cooling energy" of this break in the electron 
spectrum is (e.g. iRvbicki & Lightmanlf 19791 also, see eqs. JA4a> 
and jA5> for a derivation of this expression) 7b ~ [AoB^tjct] , 
where ijet ~ '"jot /"jet oc M is the characteristic travel time 
of the electrons. Using eq. 0, this expression becomes 7b oc 
(j}B{M, m)~^M~^. The characteristic synchrotron frequency cor- 
responding to this electron energy is 

oc B-y^ oc [MM, rh)]-^ M"" (3) 

We will assume that the spectrum has a break at Vh with pow- 
erlaw index Qx above !/b and with the standard powerlaw index 
Q = (p — l)/2 below !/b (see sketch in Fig.0. Given the expres- 
sion for Uh from eq. J3}, the flux at a fixed frequency j/x above the 
cooling break i/^ is a simple function of the flux F^q at a fixed fre- 
quency i/Q below the cooling break, the black hole mass AI, and the 
accretion rate m: 



(4) 



The uncooled, optically thin jet flux F{i'o) depends on M and m 
as follows (Heinz & Sunvaev 2003): 



F{uo) oc M^g^ {M, rn)4>c{M, rh) 



(5) 



With the expression for Fivy^) as a function of Fivo), M, and m 
from above, we have 



F{u^) oc M 



(6) 



Writing this expression in differential form, we arrive at the scaling 
relation derived in full generality in the appendix in eq. <A18> and 
eq. <AT9t 

We expect the functions 0_b(A'/, m) and (j)c{M,rn) to 



be of powerlaw type, giv en typical accretion disk scenarios; 
iHeinz & SunvaevI ^2003^ and lMerloni et all <2003h discuss the pos- 
sible functional forms of <^s and <j)c- In this case, the logarithmic 
derivatives ^xM in eq. <A18> and ^xih in eq. <A19> are constant and 
the scaling relation between Fx and M and between Fx and m as- 
sumes powerlaw form: 



Fx 



CxM = 2 — 2Qx + P + 
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Expressing the scaling relation this way is useful because p 
and Qx are actually obsen'ables. In other words, we can substitute 
observable quantities (in this case the X-ray and the particle spec- 
tral indices Qx and p respectively) for the unknown details of parti- 
cle acceleration and cooling^ . It is important to keep in mind that p 
is the electron spectral index produced by the acceleration process, 
i.e. below the cooling break and thus without any modification due 
to cooling, p is therefore not observable via the X-ray spectral in- 
dex. Rather, optically thin, uncooled radio, sub-millimeter, or IR 
emission should be used to infer p, since this radiation should not 
be affected by synchrotron losses. 

In the following we wi ll use the same canonic al functional 
forms for <f>B and <jf>c used in ( Heinz & Sunvaevl2 003). Namely, we 
assume that the underlying accretion flow is mechanically cooled"^, 
in which case the accretion flow itself is scale invariant. Then 



(10) 



This prescription is also valid for coronae above gas pr essure- 
dominated, radiatively efficient iShakura & Sunvaevlll973h accre- 
tion disks (the standard scenario for the origin of powerlaw X-ray 
emission and the launching point for jets from geometrically thin 
accretion disks, e.g. Merloni & Fabian.200Z) . We then have 

3 — Ox 



^xm 



1 — 3ax 



(11) 
(12) 



These are the expressions we will use in the remainder of this paper 

to fit the "fundamental plane" correlation. In the uncooled limit of 

Qx = (p — l)/2, they reduce to the expressions derived for the 

optically thin case from eq. jsj. In the canonical case of p = 2, we 
have Fx oc Af(3-x)/2^(5-3ax)/2^ 

Since the "fundamental plane" correlation is measured be- 
tween the radio and X-ray luminosities and the black hole mass, 
we should convert the sc aling relations to those coordinates. From 
iHeinz & SunvaevI J2003l) . the radio flux follows the scaling relation 



2p+13-(2 + p)Q 
Fr (X M 8 + 2p 



2p+13 + (p + 6)o 



(13) 



Combining this expression with eqs. jl 1> and <12> . we finally arrive 



2 iHeinz & Sunvae\< <2003h used the same approach to substitute p and the 
observable radio spectral index Or for the unknown details of the jet struc- 



ture, 
3 



m is the accretion rate at the inner edge of the disk, since this 
is where the jets are formed and is thus well defined even for 
ADIOS- ( Blandford & Beselman 1999) and CDAF- i Naiavan etal. 200(1: 
lOuataert & Gruzino\t2000ll type flows. 
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at the desired relation between Fr, M, and the synchrotron X-ray 
emission Fy^: 

Fr oc M^^^FJ"^"^ (14) 
(2p + 13 - (2 + p) Qr) (p - 1 - ax) - (2q0 



(p + 4) (2p + 1 - 3ax) 

2p + 13 + (p + 6) Qr 

(p + 4) (2p+l-3Qx) 



(15) 
(16) 



This can now be compared directly to the measured correlation co- 
efficients ^RM.FP = 0.78 and ^rx.fp = 0.6. 



3 CONSTRAINTS ON THE X-RAY EMISSION 
MECHANISM 

The "fundamental plane" correlation has the form 

log Lr = 0.78 log M + 0.6 log Lx + 7.33 (17) 

iMerloni et alJ j2003) interpreted this correlation in terms of the 
scale invariance model. From the observed correlation, it is unclear 
whether the source of the X-rays is, in fact, synchrotron emission 
from the jet, or whether it is emission from an inefficiently accret- 
ing disk; both can produce similar correlation coefficients in the 
context of the scale invariance model. The scaling relations for 
inefficient disk X-rays do pro duce a better fit to the fundamental 
plane, but lMerloni et alji2003h noted that the effects of synchrotron 
losses on the synchrotron X-ray emission from the jet had not been 
included in their treatment and might improve the fit to the funda- 
mental plane data. Given the analysis presented in the previous sec- 
tion, we can now proceed to include these effects and test whether 
they improve the fit or not. 



3.1 Spectral evidence for synchrotron cooling 

Before applying the scaling relations derived in eqs. <15> and <16t 
to the "fundamental plane" sample, we shall present a brief argu- 
ment that the inclusion of cooling into the scale invariance relations 
presented above is not just useful, but actually necessary to properly 
interpret the fundamental plane data in the context of synchrotron 
jet X-rays. 

First, we will check the global radio-X-ray spectra of individ- 
ual sources against the hypothesis that the X-rays originate in the jet 
as optically thin, uncooled synchrotron emission. In this case, the 
X-ray spectrum should continue from the X-rays to lower frequen- 
cies at the standard synchrotron spectral index Qx = (p — l)/2, 
until it reaches the turnover frequency where the jet becomes opti- 
cally thick to synchrotron self absorption (see Fig.0. At this point, 
the spectrum becomes flat. This implies that the spectrum is convex 
at every point between the radio and X-ray bands. 

Therefore, if the global spectral index Qrx is steeper than the 
fiducial uncooled X-ray spectral index ax = {p — l)/2, radiative 
cooling must have steepened the jet synchrotron spectrum at X-ray 
energies (this corresponds to point Lx3 in Fig.0. Only a cooling 
break/cutoff can reduce the X-ray flux to the observed level (even if 
the radio spectrum is optically thick)*. This constraint corresponds 
to the horizontal dashed grey line in Fig.|2| and a number of sources 
lie above this line, showing that radiative cooling must be included 

* If the X-rays originate in the disk, this limit becomes even stronger, since 
the jet-spectral index is then even steeper, thus also implying the action of 
radiative cooling. 



in the treatment of these sources. Note that relativistic beaming and 
uncorrelated radio/X-ray variability do not change this conclusion. 

On the other hand, if the global spectral index arx,FP is flat- 
ter than the fiducial uncooled synchrotron spectrum (point Lxi in 
Fig-0 all points below the dashed grey line in Fig.|2j, we cannot 
say anything about the action of radiative cooling. All we can con- 
clude in this case is that either the radio synchrotron emission must 
be self-absorbed (thus flattening the global synchrotron spectrum), 
or the X-rays must be from the disk (thus increasing the X-ray flux 
relative to the radio). 

For a sub-set of the sources in the "fundamental plane" sam- 
ple, the X-ray spectral indices ax have actually been measured, 
though for some sources with very large error bars. All of these 
sources are AGNs, as the RXTE-ASM data used in compiling the 
"fundamental plane" sample are not sufficient to determine the 
spectral slope between 2 and 10 keV. However, we note that for 
XRBs in the low-hard state, where the radio-X-ray correlation is 
observed to hold, the X-ray powerlaw spectra are typically hard - 
with 0.5 < Qx < 1.0. Most of what will be discussed regarding 
the AGN with measured Ox will also hold for these sources. 

We have plotted the global and X-ray spectral indices of all 
AGNs in the sample where ax has been measured in Fig.|2| This 
figure shows several interesting things: 



• Some sources have global spectra steeper than aix > 0.65 
(above the dotted grey line), which implies that they are affected by 
cooling or have unusually steep particle spectra. Because XRBs are 
generally less radio loud than AGNs, all XRBs fall well below this 
line (they are not shown on this plot for lack of Ox measurements), 
thus, this observational constraint on the importance of cooling is 
not available for XRBs. 

• Most of the sources in the plot have X-ray spectra steeper than 
the fiducial uncooled synchrotron spectrum of ax > 0.65 (right of 
dashed grey line). This indicates that, if the X-rays are due to syn- 
chrotron emission, (1) radiative cooling should have altered their 
spectrum or (2) the particle injection spectrum must be steeper than 
typical Fermi I shock acceleration. 

• About 25% of the sources in the plot actually fall into the area 
left of the solid black line where the X-ray spectral index is flatter 
than the global spectral index (corresponding to the dashed grey 
line through point Lx3 in Fig. 0, i.e., < arx. These sources 
have concave global spectra and the X-rays from these sources can- 
not be due to X-ray synchrotron emission from the jet, unless a par- 
ticle acceleration scheme is at work that produces concave electron 
spectra. Again, all XRBs lie well below this line, so this observa- 
tional constraint on the nature of the emission mechanism is not 
available for XRBs. We note that the error bars on Ox for all of the 
sources left of the line are large. Thus, it is prudent not to draw any 
firm conclusions from this seeming discrepancy at this point and 
await better observational determinations of Ox for these sources. 



In summary, we find that radiative cooling must be at work in 
most of the AGNs contributing to the "fundamental plane" sample, 
unless the particle spectra are intrinsically somewhat steeper than 
those produced in Fermi I acceleration and typically observed in 
optically thin radio synchrotron spectra of relativistic jets. This un- 
derlines the need for radiative cooling to be included in properly 
modeling the "fundamental plane" correlation. 
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Figure 2. Plot of the measured radio-X-ray spectral index Orx against the 
measured X-ray spectral index for the AGNs in the "fundamental plane" 
sample where this information is available, left of black line: X-ray spec- 
trum flatter than radio-to-X-ray spectrum, ax < Orx, difficult to recon- 
cile with global synchrotron spectrum; light grey area: X-ray spectral in- 
dices from best fit to "fundamental plane" correlation (see Fig.|3}; dashed 
and dotted grey lines: cooling or steep injection spectrum necessary to the 
right of and above these lines; filled dots: sources with X-ray luminosity 
Lx > 10""^ ^Eddj empty dots: sources with Lx < 10"'' iEdd- 



3.2 Fitting the scaling relations to the "fundamental plane" 
data set 

We will now compare the new scaling relations to the fundamen- 
tal plane correlation. We will assume in the following that the ra- 
dio emission is unaffected by radiative cooling. Therefore, if the 
X-rays originate in the disk, the scaling r elations from the scale 
invariance model dHeinz & Sunvaevll2003l) are unaltered and we 
will not discuss them any further, apart from reminding the reader 
that they provide the best fit to the observed correlation, provided 
that the disk emission is radiatively inefficient, (with Lx oc 
tvlerloni et al. 2003). 

If the X-rays originate in the jet, however, we must include ra- 
diative cooling, based on the arguments presented in i|3.1l The pres- 
ence of a cooling break in the spectrum was already discussed in 
321 We define the magnitude of the break as Aox = ax — (p — 1 ) /2 . 
We will assume in the following that this break lies below the X-ray 
band for all sources of interest. If this were not the case, the break 
would move through the X-ray band and produce an observable 
change in the correlation coefficients. Given the scatter in the "fun- 
damental plane" relation, we cannot rule out such a change with 
any certainty. We note that the presence of such a change would be 
an indication of jet X-rays. The dependence of the break frequency 
u\y on M, rh (or alternatively, Fx), and the radio luminosity Lr can 
be easily calculated from eqs. <A5> and <A12> and compared to fu- 
ture data. 

We will now use eqs. <15t and <16> to test whether radiative 
cooling will improve the fit of synchrotron X-rays to the "funda- 
mental plane" correlation. In the simplest possible case of contin- 
uous injection of fresh powerlaw particles with index p (e.g., by a 
stationary shock) without any adiabatic losses, the spectral break 
produced by cooling is Aox ~ 1/2. Thus, for a canonical syn- 
chrotron spectrum of index p = 2, a = 1/2, and Aox = 1/2, 
the steepened spectrum above the break has a powerlaw index 



of Qx = 1. Inserting this value into eqs. (BJ and 116} gives 
Fx oc Mm and thus ^rm ~ and ^rx = 1.42. 

This implies that cooled synchrotron emission from the jet 
follows the exact same scaling relation as efficient disk X-rays. 
As shown in Merloni et al. ( 2003), radiatively efficient disk X-rays 
are inconsistent with the observed correlation coefficients of the 
"fundamental plane". This is clear from comparing ^rm = and 
^Rx ~ 1.42 with the measured correlation coefficients of 0.78 
and 0.6 respectively. Thus, jet X-rays, in the most straightforward 
scenario of synchrotron emission (which includes cooling and the 
canonical range of particle spectral indices), are not compatible 
with the observed correlation. 

We can ask the broader question under what conditions the 
scaling relations predicted by synchrotron X-rays are actually con- 
sistent with the observed correlation coefficients. We will do this 
by fitting the "fundamental plane" data set with the correlation co- 
efficients ^RM and ^Rx from eqs. <15> and <16> . We will ignore 
upper limits in t he data and use the same merit chi-square e stima- 
tor employed in lMerloni et alj J2003h : iGebhardt et alj jioOO) . We 
also use the same error estimate, i.e., we assume isotropic errors in 
all variables and normalize them such that the reduced chisquare of 
the best fit is one. This is equivalent to equating the scatter in the 
fundamental plane to the uncertainty and is a conservative approach 
for calculating errors. 

Figure |3| shows a chi-square map for the two interesting pa- 
rameters Aax and p (marginalized over the unknown relative nor- 
malization between radio and X-ray flux). As expected, the figure 
shows that a standard cooled synchrotron spectrum with p — 2 and 
a cooling break of Aqc = 0.5 is inconsistent with the correlation 
coefficients derived from the "fundamental plane" relation by more 
than 3 sigma. The best fit actually requires a steep, unbroken pow- 
erlaw withp = 3.4l:J^;^, Aqc = Oto" and Qx = 1.2tEJj5, where 
the lower bound on Aqc 5s is set by the condition that radia- 
tive cooling cannot produce positive curvature in the spectrum over 
a large range in wavelengths. Formally, the best fit would actually 
require a negative value of Aa. 

Such a steep value of p is not in line with the particle spec- 
tra produced by Fermi I acceleration in shocks. Furthermore, this 
value of Qx is significantly larger than the mean and median of the 
measured X-ray powerlaw spectra in Fig. |2| Most of the sources 
in this figure lie to the left of the grey l-tj confidence limit on 
derived from Fig.|3| The immediate conclusion from this discrep- 
ancy is that synchrotron cooling, which should be at work in these 
sources (as argued in fj2} does not improve the fit to the "fundamen- 
tal plane" (as speculated bv ..Merloni et al.,2003.,) . To the contrary, it 
actu ally makes it worse. 

iMerloni et all |20o3) and lMaccarone et alj <2003ft discuss the 
possibility that the radio-X-ray-mass correlation of the "fundamen- 
tal plane" breaks down at high accretion rates and that different 
emission mechanisms are responsible for the X-rays above and be- 
low the break. However, repeating the analysis with a sample lim- 
ited to low X-ray luminosities (Lx < 10~''LEdd) does not change 
the conclusions reached above, only the confidence contours in 
Fig.|3lare slightly less stringent: p = S.Slfj;^, Aqx = O-l-li, 
andQx = 1.25lo;f . 

as shown by the white contours in Fig.|3] 
In fact, the X-ray spectra of the low luminosity sources (white 
dots in Fig. |2j tend to be flatter than those of the high luminos- 
ity sources (black dots). While all the observed spectral indices in 
Fig.|2|are from AGNs, the same is known to be the case in XRBs, 
which typically show harder spectra in their low luminosity state. 
Thus, most of the observed X-ray spectral indices of low luminos- 
ity sources are significantly flatter than the range in Qx allowed by 
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Figure 3. Chi-square map of the particle spectral index p and synchrotron 
cooling break Aoc , derived by fitting the "fundamental plane" radio-X-ray- 
mass correlation with cooled synchrotron einission as the source of the X- 
rays and a flat, optically thick radio synchrotron spectrum with index = 
0. Over-plotted are the the 1, 2, and 3-sigma contours (thick black lines), 
the best fit (black dot), and the canonical cooling break of Actc = 1/2 
(white dashed line). A reasonable fit requires an unbroken, i.e. uncooled 
global spectrum that is steeper than p > 3. Also shown are the same 1, 
2, and 3 sigma contours and the best fit obtained by fitting only low-X-ray 
luminosity sources with Lx ^ 10~'^LEdd (white contours and white dot). 



fitting the scaling relations for synchrotron jet X-rays to the low 
luminosity sources in the "fundamental plane" sample. 

3.3 Discussion 

This discrepancy between the correlation indices of the fundamen- 
tal plane, the predicted scaling relations, and the actually observed 
X-ray spectra suggests several possible scenarios (not necessarily 
mutually exclusive): 

3.3.1 Steep global spectra 

Fig. m indicates that the fundamental plane relation is compatible 
with X-ray synchrotron emission i/the underlying particle distri- 
bution is inherently steep, with p ~ 3.3. Such a steep, unbroken 
particle distribution can be achieved in two ways: 

a) The particle acceleration process produces inherently steep spec- 
tra withp ~ 3 (unlike Fermi I acceleration), and radiative cooling is 
unimportant even at X-ray energies. However, for radiative cooling 
to have no effect, even on the X-ray emitting electrons, the mag- 
netic fields in the inner jet would have to be very low. 

b) Particle acceleration produces spectra close to p ~ 2, but ra- 
diative cooling steepens the spectra to an index of p ~ 3 all the 
way down below radio frequencies (in which case the formalism 
above cannot distinguish between cooling and an intrinsically steep 
spectrum). However, the jet is optically thick at lower frequen- 
cies. The cooling coefficient for electrons radiating below the self- 
absorption turnover frequen cy is significantly reduced (McCrav 
Il969t iGhisellini et aljfTgsl . and an unbroken powerlaw is not a 
solution to the kinetic equation in this case. The self-absorption 
barrier can be circumvented if radiation losses are dominated by 



inverse Compton scattering. A radiative signature of this compo- 
nent must be visible, since it would dominate the radiative output 
of the jet (presumably in the X-rays). Alternatively, a somewhat 
fine-tuned adiabatic loss scheme could keep the synchrotron cool- 
ing break always below the local turnover frequency. 

However, such a steep global spectrum is inconsistent with 
the flatter radio spectral indices observed in larger scale optically 
thin jet sources (like M87, GRS 1915-1-105, 3C273, etc.) and would 
therefore require additional particle acceleration on larger spatial 
scales. It is also inconsistent wit h the standard scen ario for produc- 
ing flat radio spectra iBlandfor d & Ostrikedll97^ . And for most 
sources in Fig.|2|(and XRBs in the low-hard state), the observed 
X-ray spectral indices are much shallower than this. Ultimately, 
broad band spectra of the jet core are necessary to verify the spec- 
tral index of the optically thin part of the spectrum over a large 
frequency range and for a number of sources with different M and 
m to determine p unequivocally. 

3.3.2 Scale invariance violated 

It is possible that the X-rays of all sources in the sample are of 
synchrotron-jet origin, but that jets are inhere ntly not scale in- 
varian t and the analysis presented above and in^ein^^^^^a^ 
i2003h does not apply. However, the presence of a strong cor- 
relation and its agreement with the scale invariance predictions 
would be a coincidence. Moreover, specific models used by 
iMarko ff et al. ( 2001 ) and Markoff et al. ( 2003) (based on the origi- 
nal, scale invariant jet model by Blandford & Ko nigl.l ^ ^ , ) and bv 
Falcke & Biermann ( 1995), Markoff et al. ( 2003i). and lpalcke et alj 
( 2004) to model the correlation and the broad band spectra of indi- 
vidual sources are themselves members of the class of scale invari- 
ant models. Such a lack of scale invariance would require a more 
complex approach at understanding jet spectra and correlations. 

3.3.3 Other sources of X-ray emission 

Finally, the discrepancy can be resolved if the X-ray emission from 
at least a significant fraction of the sources in the "fundamental 
plane" sample does not originate as synchrotron emission in the 
jet, but either a) in the disk, as long as the X-ray emission from 
the disk is inefficient in the sense that oc m? , or b) as inverse 
Compton emission from the jet, with seed photons provided by the 
accre tion flow, the jet radio synchrotron emission, or the compan- 
ion (iGeorganopoulos et al. 2002: Ma,rkoff & Nowak 2004.) . 

The disk interpretation would be in line with the fact that the 
inefficient disk X-ray scenario provides the best global fit to the 
fundamental plane relation. Discussing the inverse Compton sce- 
nario is beyond the scope of this paper. It is clear, though, that the 
presence of the radio-X-ray-mass correlation will impose similar 
limits on a Compton scenario as is the case in disk X-rays (e.g., it 
might help us to determine the source of the up-scattered photons, 
i.e., disk vs. synchrotron-self-Compton). Clearly, further research 
is necessary to explore this scenario. 

Given the large scatter in the fundamental plane relation, the 
X-rays of at least a part of the sample might still be predominantly 
jet-synchrotron emission. In fact, it is quite plausible to expect that 
all three emission mechanisms will make a contribution to the X- 
ray spectrum of any given source at some level. It might well be that 
in different sources, different X-ray emission mechanisms domi- 
nate, leading to the significant observed scatter in the "fundamental 
plane". For example, the X-rays from Galactic XRBs in the low 
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hard state might be of jet syncfu'otron origin, while the X-rays from 
AGNs might originate primarily in the disk (or vice versa). In this 
case, the fact that both XRBs and AGNs fall on the same relation 
would be a coincidence^ . 

The scaling relations derived above and in Heinz & Sunvaev 
^2003^ : lFalcke & BiermannHl99.'5h can help to distinguish between 
these different types of sources in a large enough, well character- 
ized sample, which is clearly a primary objective for future research 
in this area. 



vary randomly, since the radio-X-ray correlation would otherwise 
be a coincidence. 

Thus, while we expect all of these effects to be present and to be 
important on a detailed, quantitative level, we believe that they do 
not alter the fundamental presence of the correlation and that we 
can use a statistical approach (i.e., measuring the correlation coeffi- 
cient) to learn something about the underlying nature of the emitter. 
However, all of these effects will contribute to the scatter and the 
uncertainty in any conclusion and must therefore be better charac- 
terized in future studies. 



3.4 Caveats 

Some caveats and sources of uncertainty should be kept in mind 
when interpreting these results: 

a) The "fundamental plane" sample is a relatively diverse sample, 
containing several different types of AG Ns, spanning a large r ange 
in accretion rate and black hole mass (see lMerloni et all2003l for a 
detailed discussion of the properties of the sample), mostly for the 
sake of a large number of sources and for testing the universality of 
the correlation. This is in line with the suggestion from i|3.3.3l that 
sources with different X-ray emission mechanisms may be present 
in the sample. Due to the diverse nature of the sample, some sys- 
tematic effects may be present that we have not accounted for. In 
the context of this paper, however, a reduction to the more homo- 
geneous low-luminosity sample (also used by .Falcke et aljii2004i) 
shows identical results. 

b) The X-ray spectral characteristics of XRBs change as they 
change state, and the radio and X-ray band vary on different time 
scales and with respective phase lags. This specifically affects the 
global spectral index. The same should hold for AGNs, though it 
is difficult to test due to the long time scales involved. However, 
as long as there is a simple delay between the two bands and no 
systematic trend for one band (X-ray or radio) to be in a specific 
state longer than the other, the variability should average out in a 
large enough sample and be simply another source of scatter in the 
fundamental plane. The presence of th e radio-X-ray c orrelation in 
XRBs in GX 339-4 jCorbel et alj20ol and V404 Cyg jGallo et alJ 
l2003ft indicates that variability effects do not destroy the correlation 
at least in the case of XRBs. It has to be kept in mind that the error 
bars on some of the Ox measurements from the literature shown in 
Fig.|2|are very large. Clearly, future observations with tighter error 
bars will be very helpful to strengthen the arguments made in this 
paper. 

c) The same holds for the radio spectral index: it is known that as 
a the source flux changes in XRBs, ar also varies by significant 
amounts. However, at least for AGNs we know that, on average, a 
flat spectrum is a good description of the radio core, and while the 
situation for XRBs at very low luminosities should be explored fur- 
ther, it is probably a good approximation to assume that the spectra 
are, on average, close to flat I Fender 2 001.) . Again, the radio-X- 
ray correlation indicates that the underlying spectrum is unlikely to 



4 CONCLUSIONS 

We presented an extension of the scale invariance formalism de- 
rived by Heinz & Sunvaev (2003) that takes the effects of radia- 
tive cooling into account self-consistently. The spectral index in- 
formation of the sources contributing to the "fundamental plane" 
correlation shows that the X-ray luminosity of the sources should 
be affected by cooling if it is due to synchrotron emission from 
the jet. We proved that, for the canonical synchrotron parameters 
(p = 2, Aa = 1/2), the scaling relations for X-ray and radio 
synchrot ron emission cannot r eproduce the correlation coefficients 
found bv lMerloni etai]<2003h . 

Fitting the newly derived scaling relations to the fundamental 
plane data, we showed that steep (p ~ 3.4, Ox ~ 1.2), unbroken 
spectra are required to reproduce the correlation coefficients. Such 
steep spectra are inconsistent with the observed mean and median 
spectral index of the "fundamental plane" sources with measured 
spectra. This seeming discrepancy leads us to propose three alter- 
native scenarios: a) The X-rays are of synchrotron origin and the 
entire electron spectrum has a steep index of order ~ 3, which 
would require steep injection spectra or inverse Compton cooling 
of the entire powerlaw distribution of electrons to well below ra- 
dio frequencies. However, such steep spectra are inconsistent with 
the observed range in X-ray spectral indices b) the X-rays are 
due to synchrotron emission from jets, but some of the assump- 
tions made in deriving the scale invariance predictions fail (i.e., 
some aspect of jet dynamics is not scale invariant); c) the X-rays 
from a sizeable fraction of the sources in the sample are not of jet- 
synchrotron origin and are thus most likely emitted either by the 
accretion flow itself or as Compton up-scattered radiation from the 
jet. It still plausible that a fraction of the sources (e.g., most of the 
XRBs) in the sample emit predominantly synchrotron X-rays and 
that synchrotron emission from jets contributes a small fraction to 
the X-ray spectrum of most sources. 

The analysis presented in this paper shows how a statistical 
analysis of radio and X-ray luminosities of a sample of black holes 
can be used to constrain the X-ray emission mechanism of accret- 
ing black holes and to test the scale invariance hypothesis of jets. A 
more refined sample of black holes with a better understanding of 
the sources of scatter (and the selection effects involved), and with 
tighter limits on the X-ray spectral indices is necessary to take this 
analysis to a more precise quantitative level. 



^ The fact that both populations follow very similar correlations between 
Lr and Lx for a tixed mass bin w ould be a consequence o f the fact that 
the predicted scaling relations from lHeinz & SunvaevI i20O3l) are relatively 
similar for inefficient disk X-rays and uncooled synchrotron X-rays from 
the jet. 
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APPENDIX A: FULL DERIVATION OF THE SCALING 
RELATIONS 

Al Particle transport 

In this appendix, we will derive the scaling relations for X-ray 
synchrotron radiation from scale invariant jets, including radiative 
cooling. We refer the reader to ij2|for the set of assumptions going 
into the scale invariance ansatz. Again, we assume that the particles 
are injected into the jet at some scale invariant location ro and time 
to (the latter measured in the jet frame). 

The kinetic equation describing the evolution of the particle 
distribution function / (7) as a function of time t and particle en- 
ergy 7 is (adapted from lColeman & Bicknelil988h : 



df riln(p) 
dt dt 



3 ^7 ^ ' 



2e' 



^7 



(Al) 



where an isotropic pitch angle distribution of the particles was as- 
sumed for the numerical coefficients (this assumption is not critical 
for the derivation, since the relevant numerical coefficients don't 
enter into the scaling relations). Eq. <AH can be easily expanded 
to include inverse Compton losses by adding a term proportional to 
the radiation energy density to the right hand side of eq. <A1> . 
For mathematical convenience, we substitute 



(A2) 



7 = {po/pf^^-y 

where po = p(to) is the density at the point xo = ''o/'^g where 
the particles were accelerated/injected, and t = t/to is the time 
when the particles are injected (measured in the jet frame). We also 
substitute /= (p/po)"^''V- 

As plasma travels along the jet, there is a one-to-one relation 
between its jet coordinate x and the time t, defined by 

5 = ^c/3r=^^,r(x) (A3) 
dt Tg rg 

We define xo as the acceleration location such that t{xo) ~ 1- By 
virtue of the assumption that the particle injection location is scale 
invariant, this value of xo must be independent of Al and m. We 
can then substitute the integration over t by one over x- 

The characteristic equation for synchrotron and adiabatic 
losses for each particle implicit in eq. <AH is: 



d^ 



Ao 



A ^ (—V [-^^ 

°^/3r(x) \BoJ \po J 



1/3 



-.2 

7 



(A4a) 



= ~Ao[^0r{x)r'[^B{x)f[Mx)]'^'f (A4b) 

46* 



rg4>B 



where we used eq. in the last step. 
Equation JA4at has the solution 

1 



7(x) = 



1 , 4 rx , , [V'b(x')1^K(x')]^ 
'^^ Kr(x')j 



(A4c) 



(A5) 



where 70 = 7(to)- For convenience we define 

AoT dx! [V/3r(x')l"' [^B{x')f [V'P(X')]'^' 



C(x) 



(A6a) 



= ,^^(M,m)V'5(x) 
where, by definition. 



(A6b) 



(A7) 



= 



(A8) 



and ?/>5 is independent of M and rh (but might have an implicit 
dependence on parameters such as black hole spin). ^ is simply 
the inverse of the synchrotron cutoff energy, ^ = l/7c, while ipe, 
describes how this energy cutoff in the electron spectrum depends 
on M and m. 

We can then solve eq. JA1> along the characteristics defined 
by eq. <A5> . which reduces the problem to an ODE: 

d{mx)ffmx)]] 

dx 

where d/dx is the total derivative along curves 7(x) defined by 
eq. <A5L In other words, the function 7'^/ is constant along so- 
lutions of eq. <A4at . The well known solution to eq. jAl> is then 
cKardashev 1962; Bicknell & Melrose. 1982) : 

/(7,X) = (1 -e7)-^ [Mx)\i (A9) 

for 7 5C 7c = 1/^ and the solution has a cutoff at 7c such that 
/(7 > 7c, t) = 0. For later use, we define the scaled Lorentz 
factor as the ratio of the particle energy to this cutoff: 



7 

7 = 7C = — 

7c 



(AlO) 



Using the assumption that the injected particle distribution at xo is 
a powerlaw as described in eq. J3, we can write 7(7, x) as 



/(7,x) = <t>cer'' (1 - 7)''"' [v^p(x)]^''+''^' 



(All) 



A2 Synchrotron emission from cooled particle distributions 

Following the notation in Rvbicki & LightmarJ fl979f) . we define 
the critical synchrotron frequency as 

3e 



i'c(7) 



(A12) 



16mc 

(where we have averaged over the assumed isotropic pitch an- 
gle distribution - as before, this assumption is not critical for the 
derivation below). This is the frequency where most of the en- 
ergy of the synchrotron kernel is emitted for an electron of energy 

2 

7mcC . 

The synchrotron luminosity at frequency v is then 



(A 13) 



(to) 



rf7/(7)Fs 



Ml) 



where 

Fs{x) = X I dx'K5/-j{x') 



(A14) 



is the synchrotron kernel. Substituting the expressions from 
eq. jAl H and eq. Q we get 



' xo 

»1 



2,/,(f+2)/3^P-l (A15) 



c/77-f (l-7)f-2F 
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where we have defined 

Wmccvde^ 



and 



Fo 



3e 



(A16) 



(A17) 



(where we have taken the magnetic field to be isotropically tangled, 
without loss of generality). 

Before progressing further, it is important to note several 
things about eq. <A15> 

• The low energy cutoff enters into the boundaries of the in- 
tegral over 7. However, since we are dealing with optically thin 
synchrotron radiation at high frequencies (X-rays), the low energy 
cutoff has a negligible effect on the integral, since the synchrotron 
kernel drops exponentially at low particle energies and the integral 
converges even if the lower bound is placed at zero, which we have 
done in eq. <A15> . 

• As a function of v, the integrand is entirely independent of 
M and m, which enter only into <j)B, 4>c, and 0^ in front of the 
integral. However, D does depend on M and m through (j>B and 

4>c- 

• As discussed earlier, (j>B and (j>c are simply the magnetic field 
strength B and the normalization C of the particle distribution at 
the base of the jet (both of which vary with M and m), whereas (j)^ 
describes how the cutoff energy 7c at a given position x depends on 
M and m, as defined in eq. <A7> . The fact that we can write down 
a simple dependence of 7c on M and m will allow us below to set 
the cooled synchrotron flux in relation to M and m as well. 

Following the same procedure as lHeinz & SunvaevlJiool . we 
are now in a position to derive the scaling relation between the 
synchrotron flux at X-ray energies and A/: 

dlogF^ _ d log [hp (j)B<l>cVir^) ^ dlogF^ dlogu 



aiogM 

= 2 - 2ax + p + 



aiogAf 

9 log (0c) 
dlogM 



d log v d log M 

dlog{(pB) 



+ (2p-l-3ax)- 



d log M 
(A 18) 



where we used the X-ray spectral index such that F^ cc 

_ 9 los F„ 
or Ox = — a I . 

a log 1^ 

Note that eq. (8) is fully general as long as synchrotron self- 
absorption is negligible, even in the absence of any cooling. This 
can be checked by compa ring this eq. jsj to the equivalent expres- 
sion in lHeinz & Sunvae\i' 1 2003) by substituting Qx = (p — l)/2, 
which is the appropriate expression of optically thin, uncooled syn- 
chrotron emission. 

Similarly, we can derive the scaling relation between Fi, and 
the accretion rate m: 

dlog {KP<t,B4>c€-^) 



dlogF^ 

d log rn d log rh 

d log {(t>c) 



9 log 



+ {2p-l- 3«x 



d log F„ d log V 
d log V d log m 
9 log (<^fl' 



9 log 



Cxri. (A19) 



m is simply a parameter describing the variation of magnetic field 
and particle pressure at the base of the jet at constant Af , and must 
not necessarily be directly proportional to the accretion rate. These 
two equations are the scaling laws used throughout this paper 

The superposition of cooled electron spectra from different re- 
gions in the jet might, in principle, produce any number of spec- 
tral shapes (e.g., broken powerlaws, cutoff, curved spectra, etc.). 



Eqs. {Sj and ^9), as derived in the appendix, are independent of what 
the actual X-ray spectral shape is". However, observations of syn- 
chrotron spectra from sources where we know cooling is efficient 
show either exponential-type cutoffs or, if fresh plasma is continu- 
ously supplied, broken powerlaw spectra. The latter case is appro- 
priate here, since we are observing a superposition of spectra from 
plasma that moves away from a stationary or quasi-stationary accel- 
eration region. Therefore, we can expect the cooled X-ray spectrum 
to be of powerlaw type: 



R oc M^- 



(A20) 



A3 Inverse Compton cooling 

It is easy to change the cooling mechanism in eq. <A1> to inverse 
Compton cooling off of radiation from either the central accretion 
disk or the microwave background. The latter will be entirely irrel- 
evant in the context of core dominated jets and we will ignore it in 
the following paragraph. The former can be included by substitut- 
ing Stt Umd for B± in eq. (Al\ . In this case the function (jij will 
have to be changed to reflect the dependence of the radiation energy 
density on Af and rh, which depends on the accretion scenario. In 
the case of an efficient accretion disk with (t>B = rh /Af , we have 

<j>^ oc rgM^^m (A21) 

which has exactly the same dependence on M and rh as in the 
synchrotron loss case, and the solution is unchanged. In radiatively 
inefficient flows, we have approximately 



.Af" 



(A22) 



which implies that inverse Compton is increasingly unimportant 
compared to synchrotron cooling at low accretion rates. It also im- 
plies that, while ^xM is unchanged, ^xm is changed to 

Cxm = 2p- 1/2 - Sax (A23) 
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